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ABSTRACT 

Wc present Gemini-North GMOS-IFU observations of the central starburst clumps and inner wind 
of M82, together with WIYN DensePak IFU observations of the inner 2 x 0.9 kpc of the disk. These 
cover the emission lines of Ha, [Nil], [Sii], and [Siii] at a spectral resolution of 45-80 kms^^. The 
high signal-to-noise of the data is sufficient to accurately decompose the emission line profiles into 
multiple narrow components (FWHM ^ 30-130 kms^^) superimposed on a broad (FWHM 150- 
350 kms~^) feature. This paper is the first of a series examining the optical structure of M82's disk 
and inner wind; here we focus on the ionized gaseous and stellar dynamics and present maps of the 
relevant emission line properties. 

Our observations show that ionized gas in the starburst core of M82 is dynamically complex with 
many overlapping expanding structures located at different radii. Localised line splitting of up to 
100 kms~^ in the narrow component is associated with expanding shells of compressed, cool, pho- 
toionized gas at the roots of the superwind outflow. We have been able to associate some of this 
inner-wind gas with a distinct outflow channel characterised by its dynamics and gas density pat- 
terns, and we discuss the consequences of this discovery in terms of the developing wind outflow. 

The broad optical emission line component is observed to become increasingly important moving 
outward along the outflow channel, and in general with increasing height above/below the plane. 
Following our recent work on the origins of this component, we associate it with turbulent gas in 
wind-clump interface layers and hence sites of mass loading, meaning that the turbulent mixing of 
cooler gas into the outflowing hot gas must become increasingly important with height, and provides 
powerful direct evidence for the existence of mass-loading over a large, spatially extended area reaching 
far into the inner wind. We discuss the consequences and implications of this. 

We confirm that the rotation axis of the ionized emission-line gas is offset from the stellar rotation 
axis and the photometric major axis by ~12°, not only within the nuclear regions but over the whole 
inner 2 kpc of the disk. This attests to the perturbations introduced from M82's past interactions 
within the M81 group. Finally, finding a turn-over in the stellar and ionized gas rotation curves on 
both sides of the galaxy indicates that our sight line, in places, extends at least half way through disk, 
and conflicts with the high levels of obscuration usually associated with the nuclear regions of M82. 

Subject headings: galaxies: starburst - galaxies: individual: M82 - galaxies: kinematics and dynamics 
- galaxies: ISM - ISM: kinematics and dynamics 



1. INTRODUCTION 



M82 IS tne a r cnety p e starburst galaxy 
(jO'ConneU fc Manganol [19781 : lO'ConneU et"an Tl995h 
and one of the nearest analogues to the population of 
intensely star-f orming galaxies that have been identifie d 
at high-z fe.g. iPettini et al.ll2001l : IShaplev et al.ll2003[ ). 
The current (~10 Myr) starburst activity is centred 
on the nucleus with a diameter of ^500 pc ('^30"). 
Within this region there are a number of prominent, 
high surfa ce-brightness clumps, first identified and 
labelled bv lO'Connefl fc MangaiMl (|1978[ ). From Hubble 
Space Telescope {HST) imaging, each of these clumps 
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is known to c ontain many hundred s of young massiy e 
star clusters (jO'ConneU et all Il995t iMeio et all (2005( 1 . 
and it is the energy from these clusters that drives the 
famous Ha- and X-ray-bright, kp c-scale, superwind 
(IShopbeh fc Bland-Hawtho rnI (19981 hereafter ISBH98l : 



(Stevens et al. 2003; ,Strickland fc Heckmanl(2007 

The bipolar outflow can be observed easily due to 
the galaxy's alinqst edge-on inclination (i ~ 80°; 
iLvnds fc Sai^dag^ (T963t iMcKeith etall ^M)-_The 
outflow is centred on regions A and C (ISBH98I : 
lOhvama et al.l(2"002f ). and is composed of a complex mor- 
phology of loops and filaments. Very few investigations of 
the wind dynamics have been carried out since line diag- 
nostics of the emitting plasm a are difficult to obtain with 
curre nt detec tor technology. iMcKeith et al.l ()1995( . here- 
after lMcK95( ) presented deep optical and near-infrared 
long-slit observations of emission and absorption lines 
along M82's minor axis. Clear line splitting in their emis- 
sion lines was interpreted as originating from the front 
and back walls of a cone-shaped structure. Inflection 
points in the position-velocity graph define the edge of 
a ~300 pc "energy injection region" , within which the 
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separation of the components increases rapidly (i.e. the 
outflow is accelerating). Beyond this, the optical/near- 
IR line-emitting outflow appears to re ach a ^ constant 
deprojected velocity of ~600 kms"-'^. lMcK95l inferred, 
through the difference in the strength of the red and blue 
components, that the outflow is inclined such that the 
southern cone is directed toward the observer (|McK95l . 
their fig. 5). 

What is also known is that dust in the M82 halo 
scatters emission from a central source, producing an 
extended, smooth Ha reflection n ebula (j Schmidt et al.l 
\mM iBland fc TulWl [19881 : ISBH98I ). The emission from 
this faint exponential halo is broad (^^300 kms~^), has a 
radial velocity close to ^svs, and exhibits a high [N ii]/Ha 
ratio (|Bland fc Tull v"l988'l. From their Ha hncar po- 
larisation maps. [Scarrott et al.. (199 1) inferred that the 
central illuminating source of this scattered radiation 
must be <4" (<70 pc) in diameter and located close to 
the peak of the IR emission longward of 10 fim (close 
to/within region C). This location, however, is not coin- 
cident with the dynamical centre or the position of the 
2.2 fiiB. peak usually associated with the galaxy nucleus. 
The presence of this scattered component evidently com- 
plicates any interpretation of the outflow dynamics de- 
termined from optical emission lines. 

The major drawback of long-slit studie s is the highly 
restricted nature of their spatial coverage. ISBH98I exam- 
ined the M82 wind using spatially resolved Fabry-Perot 
observations. By comparing their observations to sim- 
ple geometric models, they determined that the (south- 
ern) outflow is cylindrical out to ^500 pc, with a base 
diameter of ~350 pc, then opens out to a cone of open- 
ing angle ~25° (see their fig lib). The inner cylindri- 
cal (or pipe-like) region is associated with the aforemen- 
tioned energy injection zone. Although the line splitting 
is clear, they found that the optical emission is by no 
means distributed evenly over the surfaces of the cones, 
nor emanates smoothly from the entire starburst region. 
Clearly evident from their narrow-band im aging and line 
ratio maps (see also lOhvama et al.l l2002t ) , the southern 
wind can be traced directly back to starburst clumps A 
and C, and instead of solid cones, the optically-emitting 
gas must be excited along filam ents in the walls of the 
cone-like outflow (cf. NGC 3079: lCecil erall[200lh . Ev- 
idence that the outflow rotates in the s ame sense a s the 
disk up to heights of >1 kpc (1SBH98I : lGrevd[200l im- 
plies that rotating disk material has been entrained and 
diverted into the wind. 

Radial velocities along the major axis have been 
studied many time s in the opt ical/near-IR (e.g 
Burbidee et al.l 119641: ISaito et al " 



Sofue et al.l 



ISBH9S 



11984: McKcith et al 



WestmoQuette et al.l 


2007c; Konstantopoulos et al. 


in prep.), IR (e.g. 


Telesco et al.l 119911: iLarkin et all 


1994 lAchtermann & Lacvl 


19951) and radio (e.g. 


Shen & Lol Il995l 


Wills et al 


120001: iSeaauist & ClarkI 


20011: IWalter et al. 


120021 iRodrieuez-Rico et al.l I2004D. 



Overall, these studies have found an inner ~5" re- 
gion with a steep velocity gradient, an outer region 
>±25" of constant (or slowly increasing) velocities, and 
a transition region w here the r adial vel ocities turn over. 
iMcKeith et alTljTggsl . hereafter [McKH) were the first to 
notice that the gradient of the major axis velocity curve 
appears to increase with wavelength (from [O ii] 3729 A 



to Ncii 12.8 fim). They interpreted this as arising from 
the fact that longer wavelengths probe deeper into the 
galaxy, where for the same projected position the radial 
(line- of- sight) component of the orbital motion is greater. 

The increasing wealth of velocity meas urements have 
lead to the discovery of a ste llar bar (jTelesco et al.l 
119911: lAchtermann fc Lacvlll995t ). Through comparison 
to models, the bar was found to possess two families of 
intersecting orbits - an outer set of so-called xi-orbits 
(which form the main bar structure, and are traced 
mostly by the cooler neutral gas), and an inner set of 
perpendicularly oriented a;2-orbits contained within the 
inne r ^5'^ region, and are associated w ith the ionized 
gas (jWills et al.ll2000l: iGreve et al.ll2002D . The latter or- 
bit family arise s only in the presen ce of an inner Lind- 
blad resonance (jAthanassoulalll99^ . The interaction be- 
tween the xi and a:;2-orbit families has resulted in the 
build up of a torus of molecular and ionized gas and dust 
ident i fied in radio and IR imaging (jAchtermann fc Lacvl 
[19951: iWeifi et al.l [ 200l. Bv c ombining HST STIS spec- 
trosc opy an d imaging. IWestmo quc tte et al.l (|2007cl . here- 
after lWOTH ) suggested a possible 3D geometry for the bar 
with respect to the well-known features of the M82 star- 
burst. 

rwo73 were able to identify multiple components 
in the bright emission lines of Ha and [N ii] in 
the central regions of M82, including a very broad 
(200-300 kms~^) feature underlying the main nar- 
row component. Broad emission line wings have 
been detected in a number of other nearby starburst 
galax i es fe.g. Ilzotov et all Il996t iHomeier fc Ga ll agher 



19991 [Marlowe et all [19951: iMendez fc EstebanI il99 



Sidoh. Smith, fc CrowtheTl200fil : IVanzi et all [200fiF . 

However, due to mismatches in spectral and spatial reso- 
lution and in the specific environments observed, the na- 
ture of the energy source for these broad lines has been 
contested, resulting in a large number of possible expla- 
nations being proposed. Recent detailed IFU (integral 
field unit) studies of the ionized ISM in NGC 1569, how- 
ever, hav e shed a considerab le amount of light on this 
problem (jWestmoquette et al . 2007a^b, 2008). 

By mapping out the properties of the individual line 
compo nents (including the broad underlying compo- 
nent), [^^stmofluettiler^ ([2007a[ ) identified a number 
of correlations that allowed them to determine the likely 
origin of the broad component. They concluded that the 
evaporation and/or ablation of material from cool inter- 
stellar gas clouds caused by the impact of t he high-energy 
photo ns and fast-flowing cluster winds ([Pittard et al.l 
|2005| ) prod uces a highly turbulen t layer on the surf ace of 
the clouds ([Slavin et all 1993: Binette et al.lll999l ) from 
which the emission arises. iW07d showed that this expla- 
nation is also applicable to the broad lines in M82. They 
argued that since the high pressure ISM is highly frag- 
mented, with many small clouds well mixed in with the 
star clusters, there are many cloud surfaces with which 
the copious ionizing photons and fast winds can interact. 

Clearly, the structure and dynamics of the M82 sys- 
tem are complex. They vary both on the large and small 
scales, and are complicated by the inflows and outflows 
caused by the bar and starburst. Long slit observations 
have hitherto contributed a great deal to our knowledge 
of the rotation and outflow dynamics, but are of limited 
use in the long-term due to their restrictive linear spatial 
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coverage. We have, therefore, obtained a set of high res- 
olution, integral field, optical and near-IR spectra of the 
nuclear regions and disk of M82 at two complimentary 
spatial scales. With the Gemini GMOS-IFU we have 
observed the brightest starburst clumps and part of the 
inner wind region in detail with five 7x5 arcsec point- 
ings. With the WIYN DensePak IFU we observed four 
partially overlapping regions covering the central ~2 kpc 
of the disk. These two datasets together provide a wealth 
of information on the state of the ionized gas from both 
nebular and stellar origins, and present a unique oppor- 
tunity to study the nuclear starburst and its wider en- 
vironment. We have also obtained a set of very deep 
narrow-band images of the superwind with the WIYN 
telescope in order to investigate the large-scale structure 
of wind outflow. 

In this paper we present the IFU datasets and describe 
in detail the reduction and analysis techniques employed. 
We then focus on the gaseous and stellar kinematics (line 
widths and radial velocities), and discuss and contrast 
how the dynamics of the nuclear regions connect with 
those of the large-scale disk and the starburst system as 
a whole. In a companion paper (Westmoquette et al. in 
prep.; hereafter Paper II), we will present the second half 
of our analysis of the IFU data, focussing on the proper- 
ties of the ionized gas in the wind, including a discussion 
of the gas densities and excitations. In a forthcoming 
contribution (Gallagher et al. 2008; hereafter Paper III) 
we will present new deep narrow-band imaging of the 
wind and put forward an morphological picture of the 
whole wind system. 

In this work, we ado pt a distance to M82 of 3.6 Mpc 
(jFreedman et al.l[l994[ ). meaning 1" = 17.5 pc, a sys- 
temic velocity, Vsys, of -f200 kms~^ ()McK93l ). and 
a pho tometric major axis position angle (PA) of 65° 
(|de Vau coulc urs et aLl fTgOl: SBH98). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Gemini observations 

In February 2006 and February 2007 queue-mode 
observations using the Gemini-North Multi-Object 
Spectrograph (GMOS) In tegral Field Unit (IFU; 
lAllington-Smith et al.l[2003 ) were obtained covering six 
regions near the centre of M82 (programme IDs: GN- 
2006A-Q-38, and GN-2007A-Q-21, PI: L.J. Smith), with 
0.3-0.8 arcsec seeing (5-14 pc at the distance of M82). 
A nearby bright star was used to provide guiding and 
tip-tilt corrections using the GMOS on-instrument wave 
front sensor (OIWFS). 

Opting for two-slit mode gave us an IFU field-of-view 
(FoV) of 7 X 5 arcsccs (^^50 x 35 pc at the distance of 
M82) sampled by 1000 contiguous hexagonal fibres of 
0'.'2 diameter. An additional block of 250 fibres (cover- 
ing 2.5 X 1.7 arcsecs), offset by 1 arcmin from the ob- 
ject field, provided a dedicated sky view. We took two 
dithered exposures at each position (dither offset = 0.5 
arcsec parallel to short axis of IFU), with integration 
times between 1200 and 2400 sees each, and used the 
R831 grating to give a spectral coverage of 6100-6790 A. 
This allowed us to cover the nebular diagnostic lines of 
Ha, [Nii]AA6548,6583 and [S ii]AA6716, 6731 at a disper- 
sion of 0.34 A pix ^. Table [1] lists the coordinates, posi- 
tion angles (PAs) and exposure times for each position. 



TABLE 1 
Gemini GMOS-IFU observations. 



Pos. 


Date 




Coordinates 




PA 


Exp. Time 


No. 






(J2000) 




n 




1 


06/2/06 


og*" 


55™53?29 69° 40' 


46('8 


312 


1200 X 2 


2 


14/2/07 


OQ'' 


55™5ffi69 69° 40' 


39'.' 6 


312 


1200 X 2 


3 


14/2/07 


OQh 


55™52?31 69° 40' 


42'.' 9 


250 


2100 X 2 


4 


14/2/07 


OQ'' 


55™53?77 69° 40' 


41'.' 6 


250 


2400 X 2 


5 


15/2/07 


OQ*" 


55'"51!64 69° 40' 


34'.' 3 


250 


2400 X 2 



The GMOS spectrograph is fed by optical fibres from 
the IFU which reformats the arrangement of the spectra 
for imaging by the detector. This is composed of three 
2048 X 4068 CCD chips separated by small gaps. In order 
to remove the pixel-to-pixel sensitivity differences, and 
enable the wavelength and flux calibration of the data, 
a number of bias frames, flat-fields, twilight fiats, arc 
calibration frames, and observations of the photometric 
standard stars G191-B2B and Hiltner 600, were taken 
contemporaneously with the science fields. 

The IFU positions were chosen to cover selected areas 
of the disturbed ionized interstellar medium in the in- 
ner region of M82, near the roots of the wind outfiow; 
positions 1 and 2 cover t he bright starburst c lumps A 
and C (nomenclature from lO'Connell fc Manga no 1978), 
whereas positions 3, 4 and 5 were aligned perpendicularly 
to the major axis and cover the inter-clump region and 
the outflowing gas to the south. In Fig. [TJ we show the 
position of the IFU fields on an archive HST/ ACS F656N 
image (GO 9778, P I. L. Ho). Super star cluster (SSC) Al 
(jSmith et al.ll2006[ ) is located on the north-easternmost 
corner of position 1, as indicated with a white circle. The 
sixth IFU field was positioned over the star cluster M82- 
F located ~500 pc to the west of the nucl eus, and, since 
an an alysis of these data is presented in iBastian et al.l 
(|2007t l. will not be discussed here. 

2.1.1. Reduction 

The field-to-slit mapping for the GMOS IFU reformats 
the layout of the fibres of the sky to one long row con- 
taining blocks of 100-150 object fibres interspersed by 
blocks of 50 sky fibres. In this way, all 1250 spectra can 
be recorded on the detector simultaneously. 

Basic data reduction was performed following the stan- 
dard Gemini reduction pipeline (implemented in iraf^). 
Briefly, a trace of the position of each spectrum on 
the CCD was first produced from the fiat field. Then 
throughput correction functions and wavelength calibra- 
tion solutions were created and applied to the science 
data, before the individual spectra were extracted using 
the fiat-field trace to produce a data file containing 1250 
reduced spectra, each one pixel in width, and ordered by 
the position they were recorded on the CCD (hereafter 
the 'fibre order'). After extraction of the spectra, the x 
and y spatial units are termed 'spaxels' to differentiate 
from 'pixel', which refers to the CCD. The final steps of 
the pipeline were to clea n the cosmic-rays using LACOS- 
MIC ()van Dokkum|[200ll ) , subtract an averaged sky spec- 

^ The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by the National Optical Astronomy Observatories which is 
operated by the Association of Universities for Research in Astron- 
omy, Inc. under cooperative agreement with the National Science 
Foundation. 
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trum (computed from the separate sky field), and apply 
a flux calibration derived from observations of the stan- 
dard star (G191-B2B for the Feb 06 run, and Hiltner600 
for the Feb 07 run). Note that although the one arcmin 
separation of the GMOS sky and object fields (see above) 
was not enough to place the sky field on a completely 
dark region of the sky, the contamination from emission 
in the halo of M82 was found to be <0.5%. Separation of 
the sky spectra from the science data resulted in a data 
file formed of 1000 fully reduced spectra. 

The data were converted into standard 'cube' format 
using GFCUBE with no interpolation. To combine the 
individual dithered exposures for each position, we used 
a custom pyraf script. After the required interpolation 
was applied to combine the fields, the result was a data 
file containing 825 spaxels on which wc could begin our 
analysis. 

In order to determine an accurate measurement of the 
instrumental contribution to the line broadening, we fit- 
ted single Gaussians to three isolated arc lines on a wave- 
length calibrated arc exposure associated with position 
1, for all 1495 apertures (science-|-sky). The average 
instrumental broadening (velocity resolution; FWHM) 
of the final dataset is 1.26 ± 0.06 A, or 59 ± 3 kms^^ 
(this is co nsistent with our previous work with this in- 
strument; l^^itmofluitti^^E^llOOT^)- The wavelength 
range covered is small enough such that instrumental res- 
olution does not vary significantly from the blue to the 
red extremes. Similarly, differential atmospheric refrac- 
tion does not affect our data in any significant way. 

2.2. WIYN observations 

DensePak (jBarden et al.|[T998[ ) is a small fibre- fed in- 
tegral field array attached to the Nasmyth focus of the 
WIYN (Wisconsin, Indiana, Yale and NO AO) 3.5-m tele- 
scope. It has 91 fibres, each with a diameter of 300 fim 
or 3" on the sky; the fibre-to-fibre spacing is 400 fim 
making the overall dimensions of the array 30 x 45 arc- 
sees. Four additional fibres are offset by ~60" from the 
array centre and serve as dedicated sky fibres. The for- 
mat of the DensePak fib re array on th e sky, including 
the sky fibres is shown in lSawven ()1997f ). At the time of 
observation, there were 5 damaged and unusable fibres 
in the main array (36, 40, 46, 59 and 93) making a us- 
able total of 86. DensePak's fibre bundle is reformatted 
into a 'pseudo-slit' to feed the WIYN bench-mounted 
cchelle spectrograph. This spectrograph is part of the 
Hydra multi-object fibre positioner instrument, and uses 
a T2KC 2048 x 2048 CCD detector. 

On 14th April 2001, four fields covering the stellar disk 
of M82 were observed with DensePak. Unfortunately, 
due to a telescope malfunction, the precise coordinates 
for the fourth position were not recorded at the time 
of observation. To recover approximate coordinates we 
have matched reconstructed images in the continuum and 
Ha bands to WIYN and HST imaging. The resulting 
certainty in the coordinates for position 4 is degraded ac- 
cordingly. The 860 line mm~^ grating gave a wavelength 
range of 5820-6755 A with a dispersion of 0.46 Apix~^, 
allowing us access to a number of optical nebular lines 
(hereafter referred to as the Ha grating). The first three 
positions were also observed with a second grating setting 
with a wavelength range of 7745-9700 A and dispersion 
of 0.96 A pix , providing access to the ncar-IR Ca ii stcl- 



TABLE 2 

WIYN DensePak observations. Coordinates refer to the 

ARRAY centre (fIBRE 42). 



Position 




(J2000) 

RA 


Dec 


Exposure Time 

(s) 


Ha grating 












1 


9" 


55™ 53! 53 


69° 


40' 43'.' 14 


3 X 1200 


2 


9h 


55™ 48? 56 


69° 


40' 32'.' 20 


3 X 1200 


3 


9" 


56™ 04? 11 


69° 


41' 11'.'57 


3 X 1200 


4 


gh 


55™ 55? 7 


69° 


40' 57" 


1 X 1200 


Call grating 












1 


9h 


55™ 53? 53 


69° 


40' 43'.' 14 


2 X 1200 


2 


9" 


55™ 48? 56 


69° 


40' 40'.' 75 


1 X 1200 


3 


9" 


56™ 04? 11 


69° 


41' 11'.'57 


2 X 1200 



lar absorption lines (AAA8498, 8542, 8662) and a number 
of additional nebular emission lines (hereafter referred 
to as the Can grating). Details of the observations are 
given in Tabled A number of bias frames, flat-fields and 
arc calibration exposures were also taken together with 
the science frames. 

2.2.1. Reduction 

Basic reduction was done using the CCDPROC task 
within the NOAO IRAF package. Instrument-specific re- 
duction was then achieved using the hydra tasks, also 
within the NOAO package. The first step was to run 
APFIND on the master fiat-field exposure to automati- 
cally detect the individual spectra on the CCD frame. 
The output of this task is an aperture identification ta- 
ble containing the trace of the dispersed light from each 
fibre, which can be applied to the science frames to ex- 
tract out the individual spectra. The task dohydra was 
then run to perform the flat-fielding and wavelength cal- 
ibration. The datafile at this stage contained 90 reduced 
and wavelength calibrated spectra (one for each fibre in 
use, including the four sky fibres). 

Even with the sky fibres separated from the main ar- 
ray by 60", extended emission from the galaxy contam- 
inated, to varying degrees, each sky fibre in each of the 
positions. After some testing, we found that sky fibre 
90 of position 3 was the least contaminated. We there- 
fore formed a representative sky spectrum by averaging 
the data from this fibre from all three exposures, then 
manually removing the remaining low-level contamina- 
tion (using the 'dcblcnd' function in SPLOt). This sky 
spectrum was then subtracted from the data (after scal- 
ing to the continuum level of fibre 90 for the other po- 
sitions). C osmic-rays were cle aned from the data using 
LACOSMIC (jvan Dokkumll200ll ). before final combination 
of the individual frames using imcombine. The final 
datafiles now contained 86 (sky fibres removed) reduced 
and sky-subtracted spectra. 

In order to determine an accurate measurement of the 
instrumental contribution to the spectrum broadening 
for the two grating settings, wc selected high S/N spec- 
tral lines from wavelength calibrated arc exposures that 
were close to the Ha (Ha grating) and [Siii]A9069 (Can 
grating) line in wavelength, and sufficiently isolated to 
avoid blends. We then fitted a Gaussian to these lines 
for all 86 apertures and took the average. The instrumen- 
tal broadening is 44.7 ±3.1 km s^^ in the Ha grating and 
81.8 ±4.1 kms^^in the Can grating. 

Fig. [2] shows a reconstructed image of the three com- 
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bined DensePak fields in the line-free continuum region 
6645-6660 A, overlaid on a WIYN i?-band image of the 
nuclear regions. The inset shows the position of the 
DensePak fields on an HST+WIYN narrow-band com- 
posite (Paper III), for comparison to the Ha emission 
from the wind. Position 1 was centred on region A of 
the star burst, position 2 samples the region to the west 
of the main starburst nucleus, including the area sur- 
rounding cluster F, and position 3 was centred on region 
B. Position 4 samples the eastern part of the nuclear 
starburst region, including part of the north-south dust 
lane. In many of the fibres near the starburst nucleus, 
the signal-to-noise (S/N) is very high (>500 at the cen- 
tre of position 1). Two example spectra from each of the 
grating settings are shown in Fig. [S] 

2.3. Decomposing the line profiles 

Following the methodology employed by 
IWestmoauette et"al] ([2007a, b), we fitted multiple 
Gaussian profile models to each emission line u sing an 
IDLY^ fitting package called pan (Peak ANalvsis; iDimeol 
I2005D . to quantify the gas properties observed in each 
IFU field. A detailed description of the program and 
our fitt ing methods can be found in lWestmoquette et al.l 
(l2007aD . 

Emission lines are detected in every fibre of both the 
GMOS and DensePak datasets. The high S/N and spec- 
tral resolution of our data have allowed us to quantify 
the line profile shapes of these lines to a high degree 
of accuracy. In general, we find the emission lines to 
be composed of a bright, narrow component (hereafter 
CI; this was sometimes split into two narrow compo- 
nents, hereafter CI and C3) overlaid on a fainter, broad 
component (hereafter C2). Each line in each of the 825 
spaxels (GMOS IFUs) and 86 spaxels (DensePak IFUs) 
was therefore fitted using a single, double, and triple 
Gaussian component initial guess. Line fluxes were con- 
strained to be positive and widths to be greater than 
the instrumental contribution to guard against spurious 
results. In some regions of our DensePak data, we can 
also identify an Ha absorption component of stellar ori- 
gin (hereafter C4). Here we add an additional Gaus- 
sian absorption component, with a flux constrained to 
be negative. To fit the Ha and [Nii]A6583 lines, for each 
component we constrained the wavelength difference be- 
tween the two Gaussian models to be equal to the labora- 
tory difference, and FWHMs to equal one another. This 
was also the approach taken for fitting the [S ii] doublet. 
Multi-component fits were run several times with differ- 
ent initial guess configurations (widths and wavelengths) 
in order to account for the varied profile shapes. How- 
ever, we note that the minimisation routine employed 
by PAN is very robust with respect to the configuration 
of the fit initial guess. 

In order to determine which fits were correct, and how 
many Gaussian components best fit an observed profile, 
we applied a number of tests and filters. Firstly, to select 
which result to keep of the multi-component fits made 
with varied initial guess parameters, we simply chose the 
fit with the lowest value. Secondly, to determine how 
many Gaussian components best fit an observed profile 
(one, two or three), we used the statistical F-test. The 
F-distribution function allows one to calculate the signif- 
icance of a variance (x^) increase that is associated with 



a given confidence level, for a given number of degrees 
of freedom. The test will output the minimal increase of 
the ratio that would be required at the given confi- 
dence limit for deciding that the two fits are different. If 
the ratio is higher than this critical value, the fits are 
considered statistically distinguishable, and the one with 
a lower x^ is selected. However, this test only tells us 
which of the fits (single, double or triple component) is 
most appropriate for the corresponding line profile. Ex- 
perience has taught us that a number of physically mo- 
tivated tests are needed to, firstly filter out well-fit but 
physically improbable results, and, secondly assess which 
Gaussian profile belongs to which physical line compo- 
nent. 

For a fit to be accepted, we set the criteria that the 
measured FWHM had to be greater than the associated 
error on the FWHM result (a common symptom of a 
bad fit). Of a triple-Gaussian fit, we specified that the 
broadest component should be assigned to component 2 
(C2), and after that, the brightest to be component 1 
(CI) and the faintest narrow component to component 3 
(C3). With double-Gaussian fits, we assigned C2 to be 
the broader of the two components and CI to the other, 
regardless of their velocity difference. This consistent ap- 
proach helped limit the confusion that might arise during 
analysis of the results where discontinuous spatial regions 
might arise from incorrect component assignments. 

Experience has shown us that the errors that pan 
quotes on its fit results are an underestimate of the true 
uncertainties. This is certainly true once uncertainties on 
our post-fitting tests and filters are taken into account. 
We can therefore estimate more realistic errors through a 
visual re-inspection of the profile-|-fit after knowing which 
one was selected by our tests, and by taking into account 
the S/N of the spectrum. For the fluxes, we estimate that 
there is a 0.5-10 per cent error in CI (with the range re- 
sulting mostly from the variation in S/N levels), 8-15 
per cent on C2 and 10-80 per cent on the C3 flux. The 
presence of a third fit component may increase the errors 
in the CI and C2 fluxes by 5-10 per cent, particularly 
where C3 is faint. We estimate that errors in FWHM 
range between 0.25-2 kms'^ for CI, 4-20 kms"! for C2 
and 15-25 kms~^ for C3, whereas those for radial ve- 
locities range between 0.1-3 kms^^ for CI, 2-5 kms~^ 
for C2 and 10-30 kms~^ for C3. Similar percentage in- 
creases in errors for CI and C2 as for the fluxes may 
occur with the addition of the third component. 

The red grating of the DensePak data covered the 
Can triplet stellar absorption lines. For the fitting of 
these lines, we adopted a different method. In order to 
simultaneously fit all three lines together we employed 
the penal ised pixel fitting (pPXF ) method, described in 
detail in ICappellari fc EmsellemI (|2004f ). In brief, the 
method simply cross-correlates selected regions of a spec- 
trum with a template, resulting in a measurement of the 
best-fit radial velocity. In our case we used the NIR Ca ii 
triplet templates from lCenarro et all ()2001l ) made avail- 
able through the pPXF website^. 

Figs 13] and show a number of example Ha (-|-[Nii]) 
line profiles and best-fitting Gaussian models from a 
number of regions within our GMOS and DensePak cov- 
erage (labelled with the corresponding letters in Figs [T] 

* http:/ /www.strw.lcidcnuniv.nl/~mcappell/idl/#ppxf 
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and [21 respectively). This sample represents the variety 
of different profile shapes we find across the GMOS and 
DensePak IFU fields, and demonstrates the high quality 
of the spectra and the accuracy of the line-fitting. 

3. DYNAMICS OF THE NUCLEAR REGIONS AND INNER 

WIND 

In this section we will present and describe the 2D 
FWHM and radial velocity maps derived from the GMOS 
observations. These cover five regions near the M82 nu- 
cleus, and parts of the inner- wind as shown in Fig. [TJ 

3.1. FWHMs 

Maps of the FWHM of each of the Ha line compo- 
nents are shown in Figs. [H [7] and [H As mentioned in 
Section [2731 the Ha line profile across the five fields is, 
in general, a convolution of a bright, narrow component 
(CI) and a broad, fainter component (C2). The majority 
of position 1, for example, exhibits this kind of profile. 
In position 4 (and some of position 3), however, we see a 
split narrow component (CI and C3), and in some places 
a very broad underlying component (up to ~350 km s~ ^ ) . 
The kinematic structure of these inner regions is clearly 
complex, with evidence for increasing complexity and ir- 
regularity away from the starburst clumps, out into the 
inner wind. We will now briefly describe some of the 
details of these maps. 

In CI, the line widths are fairly uniform across the five 
fields at ~30-120 kms~^, but increase to <150 kms~^ 
in a number of localised regions. The widths of C2 
(Fig. [7]) are quite uniform across positions 1 and 3 (130- 
250 kms~^). However, there is a distinct peak in the cen- 
tre of position 2 where the line widths reach ~270 km s~ ^ . 
This is coincident with the fiux peak of region C and the 
'region C bubble' discussed below fSection 13. 2[) . In po- 
sitions 4 and 5, the FWHM variations trace north-west- 
south-east oriented patterns, and reach >300 kms^^ in 
places. The widths of C3 (Fig. [8]) are in general simi- 
lar to or less than the CI FWHMs; in the south-west of 
position 4, however, we do see the C3 line widths con- 
sistently increasing above 100 kms~^, and in position 2 
(region C) their widths are narrower than CI. Table [31 
lists the average line widths for the three line components 
over all five fields (the errors quoted represent the stan- 
dard deviation on the sample and do not refiect errors 
on the individual measurements). 

In position 4, we find a compact region of very broad- 
line gas (<350 kms~^) near the centre-north, evident in 
C2 only. Hereafter we refer to this region as the 'position 
4 knot'. An example line profile from this knot is shown 
in Fig.[4ll, showing how at this point C2 is not only very 
broad, but highly blucshiftcd (>150 kms"-'^; see below). 
We discuss its possible origins in Paper II after analysing 
the line ratio and excitation information as well as HST 
images. Wc come to the conclusion that this unresolved 
source is likely to be a star surrounded by a compact 
nebula, possibly a Luminous Blue Variable. 

In summary, wc find uniform line widths in CI and C2 
across regions A and C in positions 1, 2 and 3 (partic- 
ularly in region A). The width of the broad component 
(C2) is on average ^^200 kms^^, rising to over 300 kms^^ 
in places. In the inner wind regions it traces linear pat- 
terns parallel to minor-axis; these are particularly obvi- 
ous in C2. 



TABLE 3 

Average FWHMs (mean ± standard deviation) for the 

THREE line COMPONENTS OVER ALL SPAXELS IN ALL FIELDS FOR 

THE GMOS AND DensePak data. 



Component Average GMOS FWHM Average DensePak FWHM 

(kms-^) (kms-^) 

CI 90 ±30 116 ±41 

C2 212 ±39 224 ±52 

C3 68 ±34 103 ±40 



3.2. Radial velocities 

The aforementioned radial velocity maps are presented 
in Figs [H [TOl and [11] for the three Ha line components. 
The highest positive (redshifted) velocities in all three 
components are found in the north-east of position 1. 
However on the spaxel-scale there are distinct morpho- 
logical differences between the three component maps: 
in CI there is a smooth gradient in velocity peaking at 
~160 kms^^ and falling towards the south-west, whereas 
in C2 the velocities peak at only ~100 kms~^. The over- 
all red-blue velocity gradients seen in the three compo- 
nents over the five IFU positions reflect the large-scale 
rotation of the galaxy (including the bar in the inner 
~20"). However, overall they do not appear to be aligned 
with the photometric major axis. For example, the high- 
est negative (blucshiftcd) velocities are found in position 
5 (not 2), and again are much higher in CI ( — 150 to 
-180 kms-i) than C2 (-100 to -150 kms"!). If the 
gas follows normal disk rotation with an axis aligned with 
the major axis, then the most blueshifted velocities would 
be expected to be in position 2 since its coverage extends 
further to the west. This tilt in th e axis of rotation is in 
agreem ent with what was found bv lRodriguez-Rico et al.] 
l|2004 see their fig. 8a) from radio observations of the 
H92a recombination line, and will be discussed more in 
Section [42Tl 

Disregarding the large-scale differences due to galaxy 
rotation, the CI and C2 velocities in position 2 (region 
C), appear more uniform and smooth than in position 
1 (region A). Coincident with the peak in Ha flux and 
line widths in the centre of region C, we see a region 
of C2-emitting gas with distinctly more positive (red- 
shifted) velocities to its surroundings (by <50 kms~^). 
As alluded to above, we will refer to this as the region C 
bubble. 

The position 4 knot (identified above as a compact 
region of very broad C2-emitting gas near the centre- 
north of position 4) shows up distinctly in the C2 radial 
velocity map (Fig. [TOl) . Clearly the knot is distinct (at 
least kinematically) from its (projected) surroundings. 
As mentioned above, the nature of this knot is discussed 
fully in Paper II. 

The existence of a third Ha component provides ev- 
idence for increased dynamic complexity. In the north 
of position 1 we see CI and C3 lines of ~25-35 kms~^ 
separation. These extend into the centre of region A. 
We also identify split lines with separations from ~10- 
75 kms~^ in an area in position 3, towards the west of 
region A. The location of these split lines trace the Ha 
morphology well. The form of CI and C3 velocities in 
positions 4 and 5 (the wind positions) follow linear pat- 
terns oriented along the minor axis, as wc found for the 
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equivalent FWHM maps. In position 5, we interpret this 
splitting, in conjunction with line width (see above) and 
density (Paper II) measurements and morphologies, as 
originating from the walls of a distinct filament (chan- 
nel/chimney) in the outflow. These ideas are introduced 
in Section [5] and discussed more fully in Paper II. The 
line separations are ^25 kms~^ in the north-east, 125- 
150 kms~^ in the west and south-west, and 50 kms~^ in 
the centre-east. Throughout position 2 (region C), line 
splitting is also present. In the north-eastern half, the red 
component is stronger (hence the blue component was 
assigned to C3), whereas in the south-western half, the 
blue component is stronger (and the red component was 
assigned to C3); the line separation is ^50-70 kms^^. 

In summary, we see Ha line splitting of ~20-70 kms^^ 
throughout the starburst clumps A and C. This line split- 
ting increases with radius to over 100 kms~^ in the inner 
wind positions, and appears to follow minor axis oriented 
patterns. 

3.2.1. Major axis position-velocity graphs 

Although the type of maps described above are the 
most complete way of representing this type of 3D 
dataset, it can still be very difficult to interpret the re- 
sults in such a fashion: the data need to be 'reduced' 
further. Thus, wc have defined a number of pseudo-slits 
across our fields, from which we can extract position- 
velocity diagrams. 

The four lettered lines running north-east to south- 
west on Figs ini Uni and [TT] indicate the location of four 
0'.'4 wide major-axis pseudo-slits from which we have 
extracted the radial velocities of each Ha line compo- 
nent. Slit 'a' runs along the galaxy major axis (PA = 
65°); slits 'b', 'c' and 'd' are parallel to slit 'a' but offset 
by 2, 4 and 6 arcsecs (35, 70 and 100 pc) respectively. 
The resulting position-velocity diagrams are plotted in 
Fig.HH Ne ar-IR [S ill] and PIO radial velocity measure- 
ments from lMcK93l a re plotte d in the slit 'a' graph for 
comparison (labelled iMcUl gas'). 

A number of observations can immediately be drawn 
from these plots. Firstly, the velocities of CI (black 
points) along the major axis (slit 'a' ) follow the near- 
IR emission line measurements of lMcK93l (-1- sym- 
bols) very well. In the central, steepest part (in- 
ner 15-20"), we measure a velocity gradient in CI 
of ~15 kms^^ arcsec"^, in excellent agreem e nt with 
previo usly published measurements (|McK93l : ISBH98I : 
IW07cD. This steep grad ient is due to the bar a;2-orbits 
(|Wills et al.]|2000l : IW07cl ). The fact that we can see the 
turn-over in bar orbit velocities on both sides of the nu- 
cleus implies that we are seeing emission from the whole 
bar and thus more than half-way into the galaxy disk. 
This challenges the canonical idea of a highly obscured 
nuclear region. 

Secondly, the velocity gradient of C2 (red points) is 
much shallower than CI (black points) in all 4 pseudo- 
slits. For example, in the same region of slit 'a' as the 
CI velocity gradient was measured, we find a C2 veloc- 
ity gradient of ~9 kms~^ arcsec"^. The rotation of the 
gas emitting C2 must therefore be slower than that of the 
Cl-emitting gas within the disk and southern inner-wind. 
Furthermore, the amplitude of the CI velocity curve de- 
creases from slits 'a' to 'c' (clearly seen by comparing 
the velocities at, for example, —5" offset), although some 



contribution to this effect may result from the offset be- 
tween the major axis and the rotation axis of the gas. Slit 
'd' is dominated by small-scale velocity structure, making 
it difficult to measure meaningful gradients. However, it 
is clear that the C2 curve is still considerably flatter than 
that of CI. Wc discuss what these findings imply in Sec- 
tion o 

We now highlight some of the many small-scale dy- 
namical features seen in these plots. At —4" to —3" in 
slit 'a', we see a coherent rise in velocities traced by CI 
and C3, indicating the presence of a kinematically dis- 
tinct gas componen t. This rise, also observed in the STIS 
Ha observations of IW07d . occurs right in the core of re- 
gion A, and is visible as a red wedge-shaped region in the 
position 1 C3 velocity map fFig. [TT|) . Its existence sug- 
gests moderately fast shell expansion: assuming that the 
two components represent the two halves of a bubble- 
like structure, we can use (half) the velocity difference 
between the two components as an indicator of the ap- 
proximate expansion speed. Here we measure speeds of 
~30 kms-^ 

Large C1-C3 line splitting is seen between approxi- 
mately -1-5 to -1-15" offsets (region C) in all four slits, 
with line separations of 80-100 kms^^. Interestingly, 
here the velocity of the broad C2-emitting gas mirrors 
that of the CI gas, but always with a redshift. This large 
line splitting indicates that here we are probing a large 
column depth into the galaxy. A closed velocity ellipse 
can clearly be identified centred at +9" in slit 'b' (traced 
by CI and C3), measuring ~4" (70 pc) in extent. This 
corresponds to the region C bubble mentioned above. 
The expansion velocity of this structure is ~20 kms~^. 
Here again the C2 velocities mirror those of CI but with 
a distinct redshift. 

These small-scale features cannot simply arise from 
optical depth/obscuration effects, since wc would then 
expect the velocities to always be blueshifted relative 
to the local value (i.e. arising exclusively from the near 
side of the wind), which we do not observe. Identify- 
ing structures such as these, and the correspondence of 
the C2 velocities with these features, would be very dif- 
ficult to identify directly from the individual IFU maps, 
thus showing the advantages of plotting the data in this 
'reduced' form. 

To summarise, the velocity gradient in C2 follows a 
distinctly shallower gradient than CI - this is evident 
in all four slits (i.e. out to 100 pc from the major axis). 
The fact that we see the turn-over in the galaxy rota- 
tion curve on both sides of the nucleus, and large C1-C3 
line splitting at ^^-1-10" on the major axis, implies that 
we are seeing a considerable distance into the disk. This 
is in confiict with the high levels of obscuration usually 
associated with the nuclear regions of M82. Evidence of 
local, small-scale (~100 pc size) line splitting is found, 
implying structures expanding at 20-50 kms~^, includ- 
ing one clear coherent bubble in region C. Even though 
C2 has a different large-scale gradient to CI, it follows 
the local velocity patterns of the narrow components, in- 
dicating that the two components may be dynamically 
connected. 

3.2.2. Minor axis position-velocity graphs 

To help us interpret the evolution of the dynamics in 
the wind direction, we have also defined seven minor- 
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axis pseudo-slits; these are indicated on Figs. [H [10] and 
[TT] by the numbered lines running north-west to south- 
east. Slit 1 runs along the galaxy minor axis (PA = 
155°); slits 2-6 arc parallel to slit 1 but offset by —6, 
-3.2, +3.5, +6.0, +9.5 and +12 arcsecs (-105, -55, 
+60, +105, +165 and +210 pc) respectively. We have 
split slit 4 into 'a' and 'b' since they pass through the 
far-western edge of position 3 and the far-eastern edge 
of position 5. The widths of the pseudo-slits ('-^0'.'4) are 
indicated with dashed lines in each case; the resulting 
position-velocity diagrams are plotted in Fig. 1131 In the 
top-left panel, we have included the Ha, [Nii]A6583 and 
JS iii]A9 532 minor axis radial velocity measurements from 
lMcK95l for comparison. 

Our velocity measurements along the minor axis 
(pseudo-slit 1) only extend ~120 pc south of the nucleus. 
Here the velocities of CI and C2 are very consist ent with 
one another, and are in good agreement with the lMcK95l 
data. Between 2-6" from the major axis, we see C1-C3 
velocity separations of ~80 kms~^. With reference to the 
velocity map of Fig. [Ill we see that here the slit passes 
through one side of the region containing split lines in 
position 3. Slit 2 is positioned at the eastern-most ex- 
tent of our IFU coverage. Here the three components are 
all redshifted by <+100 kms~^ due to the rotation of the 
galaxy (indicated by the vertical dashed line) , and the lag 
of C2 behind CI is clearly reproduced along the entire 
extent of this pseudo-slit. The most prominent feature 
is the 'knot' of highly blueshifted (up to —180 kms~^) 
C2-emitting gas located at —6" (—95 pc) from the major 
axis. In slit 3, we see a gradual divergence of the CI 
and C3 velocities from —7" downwards; at —11" the two 
components are separated by > 160 kms^^. This gradual 
rise is mirrored by C2, where at —11" the C1-C2 velocity 
separation is ~120 kms~^. 

As described above, pseudo-slit 4 was split into two 
since positions 3 and 5 arc positioned close to one an- 
other but are not spatially contiguous. The discontinuity 
is represented in Fig. [T3| by splitting the plot into two. 
In slit 4a (position 3), the CI and C2 velocities arc sep- 
arated by only 20-30 kms~^ (where C2 is redder than 
CI). In position 5, however, just a few arcsecs away, we 
can identify a blueshifted and brighter narrow component 
(black CI points at —150 to —180 kms^^). We interpret 
the line emission referred to as CI in slit 4a (position 3) 
to be C3 here in slit 4b (since in assigning the line com- 
ponents we define C3 to be the faintest of any narrow 
lines identified). This highly blueshifted emission is not 
present in position 3, but identified throughout position 
5, indicating that we may miss a significant transition 
region between the two IFU fields. At —14" from the 
major axis in slit 4b, the separation of CI and C2 is 
-70 kms-i and of CI and C3 is -130 kms'^. 

In contrast, the minor axis velocity patterns in slit 5 are 
quite consistent between positions 2 and 5 (despite the 
small gap in spatial coverage). On this side of the galaxy, 
the velocity of the C2-cmitting gas again can clearly be 
seen to lag behind the Cl-emitting gas (now in the oppo- 
site sense). In position 2 (small offsets), C3 is found on 
both sides of CI; separations reach up to 40 kms~^ for 
the blueshifted C3 and 60 kms~^ for the redshifted C3 
emission. The pseudo-slit was purposefully positioned to 
pass through a part of IFU position 5 that contains C3 
emission, and we find C3 to be redshifted with respect to 



CI by <70 kms^^. The spatial coverage of slit 6 is fairly 
restricted; this pseudo-slit is located at the western-most 
edge of our IFU coverage, and as such, all the velocities 
are blueshifted by ~— 100 kms""'^ due to galaxy rota- 
tion (again highlighted by the vertical dashed line) . The 
presence of high velocity C3 demonstrates the presence 
of disturbed kinematics in this part of the disk. 

As mentioned above, our spatial coverage of the minor 
axis is very limited; we will therefore make use of the 
lMcK95l minor axis data in order to comp are the r esults 
from the different slits. At >-9" offset, lMcK95l iden- 
tify two rapidly diverging line components, and interpret 
them as being emitted from the front and back walls of 
the large-scale wind 'cone', where the increasing velocity 
separation is taken to be evidence of a rapidly acceler- 
ating outflow. This diverging pattern is clearly mirrored 
in slit 3 (only ±~70 pc from the minor axis), where we 
can equate their "red and blue" line components to our 
CI and C3. In this slit, however, the radius at which CI 
and C3 begins to separate appears to be much smaller 
(<7"). In slits 2 and 5, the situation is reversed, and we 
see C1-C3 line splitting at small radii (of ~40 kms^^). 
These results attest to the existence of a complex kine- 
matic environment here at the base of the large-scale out- 
flow cone, with many overlapping expanding structures 
located at different positions. 

Where we observe narrow-line splitting, we find that, 
in general, the strongest component (CI) on the eastern 
side is the redder one, whereas on the other side of the 
disk, the strongest component is the bluer one. Thus 
in the outflow cone interpretation, the back wall of the 
cone is brightest on the eastern side (which is moving 
away from us due to the galaxy rotation), whereas the 
front wall is brightest on the western side. Furthermore, 
all six minor axis pseudo-slits show evidence for increas- 
ingly blueshifted line emission with distance south of the 
major axis. This is consistent with previous observa- 
tions of the minor axis velocities (lMcK95l : [SBH98l : iGrevd 
|2004| ) , and is due to the inclination angle of the galactic 
disk (oriented such that the southern outflow is directed 
toward the observer). 

3.2.3. Summary: GMOS radial velocities 

There is a great deal of information contained in the 
GMOS radial velocity maps for the five IFU positions 
and three line components. Here we briefly summarise 
our main findings from these data. Firstly, we find that 
the (major axis) velocity gradient in C2 follows a dis- 
tinctly shallower gradient than CI, and this is evident at 
distances of up to 100 pc from the major axis. However, 
even though C2 exhibits this large-scale gradient differ- 
ence, it clearly follows the local velocity patterns of the 
narrow components, indicating that the two components 
are somehow dynamically connected. We also find that 
the amplitude of the CI rotation curve decreases with 
distance from the disk midplane. 

Evidence of Ha line splitting (i.e. presence of C3) is 
found throughout the five regions observed. In clumps A 
and C, its amplitude is ~20-70 kms~^; this rapidly in- 
creases to >100 kms^^ in the inner wind, where the mor- 
phology of the split line regions follows radia lly oriented 
patterns. In t he outflow cone interpretation (jGotz et al.l 
[l990tlMcK95h . this increase reflects an accelerating out- 
flow, where the line emission originates from the front 
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and back walls of the flow. Local, small-scale (~100 
pc size) coherent features in the CI and C3 lines are 
interpreted as evidence of structures expanding at 20- 
50 kms^^ within this outflow. Now that we can identify 
and map each individual line component, it is clear that 
this central region, at the base of the large-scale out- 
flow cone, represents a chaotic, complex kinematic en- 
vironment with many overlapping expanding structures 
located at different radii. The interpretation of the M82 
wind originating from a singular, organised "outflow cone 
base" appears to be a signiflcant over-simplification. 

3.3. FWHM vs. radial velocity 

An additional way of 'reducing' the data further to aid 
interpretation of the results is to plot the Ha FWHM vs. 
radial velocity for each IFU field, as shown in Fig. [HI 

Plotting the data in this fashion clearly shows the dis- 
tributions in width- and velocity-space for each of the 
components in each of the fields, hence illustrating the 
distinction between CI (and C3) and C2. Within each 
field there is little or no systematic offset in velocity be- 
tween the narrow (CI and C3) and broad (C2) lines. 
This is very similar to wh at we found in NGC 1569 
(|Westmoquette et al.l [2007al lbl) . The exception is posi- 
tion 4, where the CI lines are consistently redshifted com- 
pared to both the C2 and C3 points. This position also 
contains the position 4 knot mentioned above. Points as- 
sociated with this knot arc enclosed by a dashed ellipse 
in Fig. [H 

In positions 1 and 3, C3 lies very much in the same ve- 
locity range as CI. In position 2. however, there appears 
to be two populations of C3 lines, one with similar radial 
velocities to CI and one that lies to the redward side. 
In positions 4 and 5, the population of C3 points lie dis- 
tinctly to the blueward and redward sides, respectively, 
of the bulk of the CI lines. This means that on the east- 
ern side (position 4) the redshifted narrow component 
is stronger (hence assigned to CI), but on the western 
side (position 5) the blueshifted component is stronger, 
reiterating what we found above. 

4. DYNAMICS OF THE DISK 

In this section we present and describe the 2D FWHM 
and radial velocity maps derived from our DensePak ob- 
servations. In contrast to the small-scale, high-resolution 
GMOS-IFU observations that focus on just the nuclear 
regions, these data cover most of the inner ^2 kpc of the 
M82 disk out to disk heights of >500 pc (as shown in 
Fig. [2]), albeit at considerably lower spatial resolution. 
The spectral resolution and S/N of the two datascts, 
however, are comparable, thus allowing us to make mean- 
ingful comparisons between the small- and large-scales as 
sampled by the two spatially overlapping datascts. 

As described in Section 12.31 the shapes of the emis- 
sion line profiles throughout much of the four DensePak 
fields are a convolution of a bright, narrow component 
(CI) and a broad, fainter component (C2). In some re- 
gions, however, we can identify two narrow components, 
often both superimposed on the broad component. We 
refer to these components as CI and C3, and define C3 
as the fainter of the two regardless of their wavelength 
relationship fSection l2.3|) . In some places we can identify 
an Ha absorption component (of stellar origin), and fit 
it using a single Gaussian component (see, e.g., Fig.[5]i). 



To measure the velocity of the Can triplet absorption 
lines we used the pPXF method (again see Section [2T3|) . 

4.1. FWHMs 

Maps of the Gaussian FWHM (corrected for instru- 
mental broadening) of the Ha CI and C2 emission com- 
ponents, [Sii]A6717,6731 and [Siii]A9531 are shown in 
Fig. [ini The first thing to note is that the Ha and [S ii] 
maps are very similar, thus providing a certain level of 
confidence in our results. The FWHM of CI of Ha and 
[Sii]^ ranges between 60-170 kms~^. The widths of Ha 
C3 (not shown) also fall in this range. In many cases, the 
spaxels exhibiting broader than average Ha (or [Sil]) CI 
widths are also those that do not have a correspond- 
ing C2 detection. This suggests that in these spaxels 
C2 does in fact exist, but at too low a level to have 
been confidently fit. and its presence biases the CI fits 
resulting in artificially elevated FWHMs. Where it is de- 
tected. Ha (and [Sii]) C2 has widths in the range ~200- 
300 kms~^, with the narrowest lines found near the ma- 
jor axis. Where we have common spatial coverage with 
our CMOS data, our width measurements are in very 
good agreement. At distances of 400-500 pc south of re- 
gions A and C (in positions 1 and 2) , the FWHM of Ha 
and [Sii] C2 reaches <370 kms^^. This region is consis- 
tent with the location of the southern inflection p oint in 
the w ind outfiow velocities identified by .McKcit h et al.l 
(|l995f l. These infection points are indicated on Fig. [15] 
with short thick lines. A similar increase in widths in 
the north is not convincingly seen, although taking the 
Ha and [Sii] maps together, the line widths certainly do 
increase away from the disk towards the north. 

The [Siii]A9531 FWHM maps (bottom two panels in 
Fig.fTS]) arc distinctly different to the Ha and [S ii] maps. 
Overall [Sili] CI is signiflcantly broader than the cor- 
responding Ha or [Sii] CI. With widths in the range 
100-250 kms~^ and the narrowest lines found along the 
major axis; the [S ill] CI map in fact has more in common 
with the Ha and [Sii] C2 maps. In Fig. [THlwe compare 
the FWHMs of all the emission components identified in 
Ha, [Sii], [Siii], and P9 for IFU positions 1 and 2. In 
this representation, the broader widths of the [Siii] CI 
lines compared to the equivalent Ha (and [Sii]) CI lines 
is obvious (in fact the mean Ha CI FWHM in positions 
1 and 2 is 99 kms^^, whereas the equivalent [Siii] CI 
average is 127 kms~^). 

The elevated [Siii] CI widths could result in three 
ways: (1) [S ill] is probing deeper within the galaxy where 
the Cl-emitting gas is more turbulent ([Siii] is emitted 
in the near-IR so is less affected by dust extinction); (2) 
[S III] is probing a hotter gas phase (S^"*" has a higher 
excitation energy than Ha or S"*"); or (3) the presence of 
an imidentificd broad component has biased the [S ill] CI 
fits to artificially high widths (in a similar way to that 
described above) . To investigate which is the most likely 
cause, we have performed a number of tests. If the in- 
creased widths of the [S ill] line arc due to a depth effect 
(i.e. explanation 1), then the near-IR Paschen lines would 
also exhibit increased widths compared to Ha. Fig. [16] 
shows that this is not the case, and the mean P9 FWHM 

9 Recall that the [S ll] A6717,6731 doublet were each fit with mul- 
tiple Gaussian components where the widths of each line were con- 
strained such that FWHM([Sn]A6717) = FWHM([S ii]A6731) 
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over positions 1 and 2 is 102 kms^ , in good agreement 
with that of Ha. Furthermore, the major axis radial 
velocity gradients of the [Siil] CI and Ha CI lines are 
very similar, supporting the conclusion that they both 
originate at approximately the same depth within the 
disk. Could the broadening be due to biassed fits? In 
Fig. [T7] we show two example [Siii]A9531 line profiles 
together with the Gaussian models needed to fit their 
shapes, one originating far from the disk and exhibiting 
two separated components (left panel), and one high S/N 
line from near the major axis (right panel). Comparing 
these to the example line profiles shown in Fig. SI it is 
clear that the near-IR [S ill] lines do not share the same 
'bright narrow and broad underlying' shapes as the opti- 
cal emission lines. Even at high S/N, there is no evidence 
for broad underlying emission, and the profile is well-fit 
by a single Gaussian. It appears, then, that the [S ill] 
lines are intrinsically broader than their Ha equivalents, 
and there is no undetected broad emission affecting our 
fits. We can thus rule out explanation 3. This leaves 
us with explanation 2, that [Siii] is probing a hotter, 
more turbulent gas phase than Ha CI. We discuss the 
consequences of this is Section [5l 

4.1.1. Summary: DensePak FWHM results 

In Ha and [S il] , the narrowest C2 lines are found near 
the disk midplane, with the line widths increasing with 
height above the disk. The broadest C2 widths a re found 
near southern outfiow inflection point found by lMcK95l . 
but there is not an equivalent in the north. The [S ill] CI 
widths are much larger than the equivalent Ha CI widths 
(100-250 kms~^ compared to 60-170 kms^^), and the 
distribution of [S ill] CI widths closely matches that of 
Ha C2 (i.e. narrowest along major axis). After testing 
a number of possible reasons for this, we conclude that 
the [S III] emission arises from a hotter, more turbulent 
gas phase than that responsible for the Ha CI emission, 
but, unlike Ha C2, originates from the same depth along 
our linc-of-sight through the galaxy as Ha CI. 

4.2. Radial velocities 
4.2.1. Ha radial velocities 

Radial velocity maps for the four Ha line components 
(three emission and one absorption) arc shown in Fig. 1181 
The large-scale velocity gradient due to galaxy rotation is 
immediately obvious in the CI and C2 maps. The highest 
positive (redshifted) velocities arc found in CI in the east 
and north-east of position 4 (~150 kms^^), and highest 
blueshifted velocities of ~ — 120 kms"-'^ arc found in the 
south-west of position 1, extending into the centre and 
south of position 2. C2, where detected, closely follows 
the pattern of CI velocities, but with a much shallower 
velocity gradient. This mirrors exactly what we found 
with CMOS for the nuclear regions, meaning that the C2- 
emitting gas in the entire disk is rotating at a slower rate 
than CI. The significance of these findings is discussed 
in Section [5l 

Although the C3 velocity map is fairly noisy, we can 
still discern the effects of disk rotation. The presence of a 
kincmatically distinct C3 component can be interpreted 
as evidence for disturbed kinematics and/or expanding 
gas motions. The most redshifted C3 lines are found in 
position 3, in the centre of region B; here it is separated 



from the CI fine by 50-80 kms^i. C1-C3 velocity dif- 
ferences of >50 kms""'^ can also be seen in many other 
regions across the disk. 

All three line component maps show again that the axis 
of rotation of the Ha-cmitting gas is not aligned with 
the photometric major axis. The offset is such that the 
most rapidly receding (redshifted) gas is located to the 
north of the major axis, and the most rapidly approach- 
ing (blueshifted) gas lies to the south. As mentioned 
above, a similar rotatio n axis tilt was found i n the cen- 
tral regions of M82 by iRodriguez-Rico et al.l (|2004 see 
their fig. 8a) from observations of the H92a recombina- 
tion line. We believe that this offset has not previously 
been noted in the optical, and its detection re-emphasises 
the advantages of integral field spectroscopy. 

4.2.2. Near-IR emission and absorption line radial velocities 

Fig. [TOl shows radial velocity maps in the [Sili]A9531 
emission line (derived from simple Gaussian fits) and the 
Can triplet absorption lines (derived from the pPXF 
method) of stellar origin. 

The [S III] velocities mirror the Ha velocities very well, 
as expected. Like Ha, the major axis of rotation ap- 
pears to be offset with respect to the photometric ma- 
jor axis. In contrast, the Can map (lower-right panel) 
shows that the stellar rotation axis is consistent with 
the photometric major axis. We measure (by eye) the 
PA of the gaseous rotation axis to be PA = 53° (in good 
agree ment with the observations of lRodriguez-Rico et al.l 
|2004[) . This is indicated as a dashed line on the CI map 
of Fig. [THl The offset, therefore, between the photomet- 
ric/stellar rotation axis and the gaseous rotation axis is 
~12°. Interestingly the molecular CO gas does not ex- 
hibit this rotation axis tilt within the central regions, and 
instead appears to rotate with the same axis as th e stars 
(ISofue et al.lll992t ISeaquist fc Clark, 2001. : Walter et all 

4.2.3. Major axis position-velocity graph 

In order to aid interpretation of these data, we have 
extracted velocities for each component (emission or ab- 
sorption) of Ha, [Siii]A9531, P9 A9229 and Can, within 
a pseudo-slit 7" in width, along the galaxy major-axis. 
These are shown in Fig. [20l where we have split the 
graph according to the gaseous or stellar origin of the 
lines. For reference, we have also pl otted the [NiiJA6583 
major ax is velocit y measurements of lCastles et al.l ()1991l 
here after |CMG9 iD and the [Siii], P8, and Can velocities 
from|McK9i 

The spatial-resolution of the DensePak IFU is much 
coar ser than previously published long-slit measurem ents 
(e.g. ICMG91I . lMcK93l . IWestmoauette eiTall [2007a) and 
the CMOS data presented in Section [31 We will there- 
fore not discuss these results with respect to the nuclear 
regions (inner ±10"), besides noting that they are in 
very good agreement. However, the advantage of a much 
higher spectral resolution and good S/N has meant that 
we have been able to identify and measure multiple indi- 
vidual line components. We can therefore begin to disen- 
tangle some of the complexities of the position-velocity 
diagram. 

A great deal of local variation in the gas velocities is 
seen, particularly outside the X2-orbit region. For exam- 
ple, at '^-f 20" there is a dramatic fall in radial veloci- 
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ties, and, with the inclusion of the ICMG91I and lMcK95l 
data, we can identify both small- and large-scale "wig- 
gles" between —20" and —80". As mentioned above, the 
shallower velocity gradient of Ha C2 (red squares) com- 
pared to CI (black circles) is clearly evident out to ±50". 
Although C1-C3 line splitting is seen within the central 
±30", the most striking result is the 50-80 kms~^ split- 
ting seen at distances of more than —50" from the nucleus 
(i.e. in region B). Clearly this intermediate age region 
(few 100 Myr; ISmith et alll2007l : iKonstantopoulos et all 
|2008() . which is not part of the current starburst, is still 
producing high velocity gas. 

Our Ha stellar absorption velocities agree very well 
with the Call ste llar absorption line me asurements pre- 
sented bvlMcK93land lGreve et"al] (l2002f) a nd the CO ve- 
locity measurements of lSofue et al.l ~( 1998( 1. At the radii 
at which we can identify the Ha absorption component 
(C4), the corresponding emission-Hne (Ha, [Siii], P9) 
velocities are l ower, but again i n good agreement with 
the hterature (|CMG91l : [McK93f l. The extreme red and 
blueshifted C an velocities at dz5" (as sociated with the 
bar x^-orbits: IWills et al.|[2000l : IW07d ) are not observed 
because our spatial resolution is too coarse. At radii 
>±10", outside of the central a;2-orbit region, the radial 
velocities of Can exhibit a consistent, but gradual in- 
crease (in the blueshifted sense). The equivalent increase 
on the other (eastern) side of the galaxy is not covered by 
our DensePak pointings. At radii greater than —50" (for 
which, again, we do not have equivalent spatial coverage 
on the opposite side), the Can velocities remain very 
consis tent at ^110 kms^^. Evidence from iGreve et al] 
()2002f ) suggests that beyond this radius, the Can veloc- 
ities remain constant or increase slightly. 

The significant differences in stellar and gas radial ve- 
locities at radii <— 30" and >+50", together with the 
fact that the gas rotation curve appears to fall back to 
zero velocity at +60/-75 arcsecs, can now be understood 
as resulting from the offset between the gas rotation axis 
and the photometric and stellar rotation axes. The peak 
redshifted and blueshifted velocities in the emission line 
and stellar velocities are numerically very similar, but 
since they are spatially offset, our pseudo-slit, oriented 
according to the photometric axis PA, misses the gaseous 
radial velocity peaks. 

4.2.4. Minor axis position-velocity graph 

Fig. [21] shows the minor axis position-velocity graph 
for the emission line velocities only (split into the Ha 
and near-IR measurements), togeth er with the Ha, 
[Nn]A6583 and [Sni]A9532 data from|McK9i The in- 
flection points in the minor-axis position-velocity plot 
found bv lMcK95l lie at ~±20" (±350 pc). Rccah that, 
in their interpretation, the region within these markers 
is the "energy injection zone" where the outflow is accel- 
erated to its terminal velocity. 

Although ou r measur ements arc broadly in agreement 
with those of lMcK95l . there are a number of differ- 
ences. For example, in the north we do not see the 
100-150 kms^i blueshifted branch or the 50-100 kms^^ 
redshifted branch (although there are hints of both in 
C3). Here CI of the Ha line does become increasingly 
blueshifted, but only up to ~60 kms"^ (at +23"). The 
Ha C2 velocities diverge consistently from CI from '-^12", 
becoming increasingly redder. [S ni] and P9 remain be- 



tween and —50 kms~^ out to >±20". Velocities of 
Ha CI and C2 in the central disk region (within ±5") 
are fairly consistent with one ano ther, but are, on aver- 
age, redder than those of lMcK9"5l by ~20 kms~^. Here 
we identify Ha line splitting (C1-C3) of 40-50 kms~^. 
In the sout h Ha CI broadly follows the red branch of 
the lMcK95l data, together with [Sni] C2. In this region, 
[Sni] CI is redshifted by 100-150 kms~^, and appears to 
be more consistent with the blue branch. Throughout the 
entire southern region. Ha C2 is consistently blueshifted 
with respect to CI. The differences between our results 
and that of lMcK95l may arise from a number of rea- 
sons, including slight pointing or seeing differences, or 
our more detailed line profile fitting, including the iden- 
tification of the broad underlying component, C2, which 
reveals a more complex situation than could be seen in 
previous data. 

4.2.5. Summary: DensePak radial velocities 

To summarise this section, firstly we confirm that the 
axis of gaseous rotation is offset from stellar rotation axis 
and photometric major axis, and extend this finding to 
the entire inner 2 kpc of the disk. We measure the offset 
to be ~12°. Furthermore, we confirm that the major 
axis velocity gradient of the Ha broad component, C2, 
is much shallower than that of CI, and that this is true 
out to ±50" (±900 pc). 

Line splitting is observed in the narrow component of 
Ha and [Sni] of 40-50 kms""'^ throughout the central 
disk region. Strikingly, we also find splitting in Ha of 50- 
80 kms~^ in the centre of region B. Our velocities along 
th e minor axis are in broad agreement with those found 
bv lMcK95l . In the northern outflow. Ha C2 is redshifted 
with respect to CI, whereas in the south it is blueshifted. 
The C1-C2 differences range up to 80-90 kms~^. 

5. DISCUSSION 
5.1. The optical emission line components 

As set out in the introduction, we have recently begun 
to build up a physical model that attempts to explain 
the origins of the different optical emission line compo- 
nents in winds driven by starbursts based on our obser- 
vation s of the NGC 15 69 star burst (jWestmoquette et al.l 
I2007airbl . [20081 see also IWOTd ). The data presented here 
for M82 fits well with our ideas, and is of sufficient qual- 
ity to allow us to develop this basic model further. 

As a reminder to the reader, our model accounts for the 
narrow components of Ha (and the other low-ionization 
optical emission lines), referred to as CI and C3 (see 
Paper II), as originating from the photoionized gas in 
the disk. This gas is fairly turbulent (typical FWHM 
50-120 kms~^), and shows copious evidence for localised 
expanding structures (shells/bubbles) through its spht 
lines (i.e. where C3 is observed) and the morphology seen 
in HST images. In contrast the broad component (C2) 
represents emission from highly turbulent mixing layers 
on the surface of denser gas clouds, set up by the impact 
of high-energy photons and fast-flowing winds from the 
massive star clusters. Since material is easily evaporated 
and/or ablated from these turbulent layers, C2 should 
trace locations of mass-loading sites within the wind. 

M82 is a dynamically complex system, perturbed sig- 
nificantly by its gravitational encounter with M81 some 
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2 X 10^ yrs ago (|YurJll999f ). This complexity, however, 
is not easily disentangled due to the galaxy's high incli- 
nation, strong and variable obscuration, the presence of 
a bar, and a powerful starburst driving a large-scale out- 
flow. Each sight line towards the centre of the galaxy will 
traverse considerable material out of the galaxy's mid- 
plane before intersecting and possibly passing through 
a disk/bar star- forming region. Furthermore, both radio 
recombination line studies and our optical measurements 
show that the ionized gas and stellar rotation axes are 
offset by ^12°. The impact of these geometrical factors, 
in combination with substantial local structure (e.g. the 
pseudo-slit 'b' in Fig. [T2|) . naturally complicate the in- 
terpretation of the velocities. 

5.1.1. The narrow components (CI and C3) 

The generally good agreement between stellar veloci- 
ties and Ha CI along pseudo-slit 'a' (Fig. [T2|) for clump 
A, implies that CI originates in the disk of M82. This 
is also consistent with the presence of compact, bright 
emission structures in direct narrow band images of 
clump A that are suggestive of photoionized surfaces 
in or near the young star clusters (see Fig. [1]). Addi- 
tionally, the fact that the widths of Ha CI and C2 re- 
main fairly constant throughout the cores of the star- 
burst clumps is consistent with the known uniformity of 
the ga s excitation and ionization pa rameter in these re- 
gions (jForster Schreiber et al.l[2003l . Paper II), and with 
the idea of a highly fragmented and well-mixed ISM 
capable of rapidly adjustin g to environm ental changes 
(|Lord et al.. 1996.: ' Fors ter Schreiber et al.i 2003: W07c). 

The detection of C3 indicates strong line-splitting in 
the narrow component arising from dynamically expand- 
ing gaseous structures. Neither finding line splitting 
throughout the starburst core of M82, nor the compara- 
tively low implied expansion velocities, are unsurprising 
since this is a very dense but dynamically active region at 
the roots of the large-scale superwind. In addition since 
many of these features have small sizes (see Section 15. 2p 
they will often crowd along sight lines giving rise to the 
chaotic appearance and over- turbulent velocities of the 
M82 disk. 

Our [S III] line measurements present additional con- 
straints on the emission line models. The differences 
in dynamics between this higher excitation, near-IR line 
and the optical emission lines (Ha, [Sii], etc.) indicate 
that [S III] is probing a hotter (and therefore more tur- 
bulent) gas phase, co-located with Ha CI. This suggests 
that we may be seeing three distinct gas phases: one 
originating from the relatively quiescent ionized gas in 
or near the star clusters (Ha - and [Sii], [Nil], etc. ~ CI 
and C3); one originating from the same place but from 
a slightly hotter gas phase ([Siii] CI); and one from the 
turbulent mixing layers on the surfaces of ISM clumps lo- 
cated in moderate-to-low density regions exterior to the 
starburst core in the inner wind (Ha C2). 

5.1.2. The broad component (C2) 

One of the most intriguing results of this study is the 
difference in disk-wide dynamics found between the nar- 
row and broad ionized gas components. Seen in both our 
high-resolution CMOS observations, and throughout the 
inner disk with our DensePak observations. Ha C2, in 
general, exhibits a shallower major axis velocity gradient 



than CI. Investigating why this is the case is one key to 
understanding where C2 is emitted from and hence what 
it represents physically. 

Along each sight-line through the galaxy there 
is a combination of emission due to reflection off 
dust, as indicated by po l arisation measuremen t s (e.g. 
Notni fc Bronkallal 119841: IScarrott et all Il991l : iJonesI 
20001) ■ as we l l as sources that are seen directly. 
Bland fc TuUvl ()1988[ ) discovered a weak, uniform, broad 
(300 kms~^) emission component in the M82 halo which 
they associated with scattered light from the nuclear re- 
gion. Our high resolution, spatially resolved observations 
indicate that the broad C2 emission may in fact be com- 
posed of a local and a scattered light component. How- 
ever, we believe that the scattered light component must 
be a minor contributor within the regions we have ob- 
served (the starburst core and the inner disk plane) for 
the following reasons: the C2 velocities vary locally and 
are imprinted with the large-scale galaxy rotation pat- 
tern, and the line widths, densities and line ratios vary 
locally and with radius (Sections [3] and 21 Paper II). 

A further clue to help us understand the relationship 
between CI and C2 comes in the form of the filament 
described in Section [XH 

5.2. Observations of an individual outflow channel 

Position 5 of the CMOS observations covers a partic- 
ularly interesting region of the inner outflow - what ap- 
pears to be a distinct outflow channel directed away from 
the starburst clump. By examining the gas properties 
here, we can begin to better understand the relationship 
between the (Ha) narrow and broad components. Along 
the length of this filament, the narrow component is split 
with a velocity difference of up to 60 kms^^. Lateral cuts 
across the filament reveal velocity ellipses traced by CI 
and C3, whereas C2 appears at rest with respect to the 
centre of C1-C3 expansion. The edges of the filament are 
defined by broader C2 widths, higher C2 densities, and 
low CI densities ('^700 cm~^ and '^200 cm""^, respec- 
tively; Paper II), whereas in the centre of the filament, 
CI exhibits a higher density (~700 cm"'', particularly at 
the end nearest the clumps; Paper II) and the width of 
C2 decreases. These facts suggest that this filament is 
a distinct outflow channel through which hot gas is es- 
caping, interacting with the cooler ISM, and entraining 
it into the flow. This structure is direct evidence for an 
interaction between the hot and cool gas phases. 

This filam ent bears a remarkab l e simi larity to to those 
predicted bv lTenorio-Tagle et al.l (|2003l ). who attempted 
to model the outflow resulting from multiple, closely 
spaced super star clusters (SSCs). They found that the 
interaction of winds from the neighbouring SSCs leads 
to a network of filaments that originate near the base of 
the outflow. These structures represent the condensed 
gas formed through a combination of stationary oblique 
shocks and crossing shocks that occur between the neigh- 
bouring outflows, and may suffer strong radiative cooling 
depending on the local values of density, temperature and 
metallicity (and thus emit strongly at both optical and 
X-ray wavelengths). 

The scale of the M82 starburst and the sheer number 
of SSCs that it contains raises the question "why do we 
not see clear evidence of outflow ch annels everywhere?" . 
The answer may lie in the fact that iTenorio-Tagle et al.l 
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model was set up with favom'able conditions for the pro- 
duction of coherent channels (the SSCs all had same age 
and mass, and were distributed at set separations on sin- 
gle plane). The reality of irregular spaced clusters of dif- 
fering masses and ages is likely to produce a much more 
chaotic and diffuse flow, although channels similar to the 
ones predicted clearly have formed. 

The high CI densities near the base of the channel im- 
ply that this is where coUimation occurs, and the expan- 
sion of the narrow component-emitting "walls" results in 
the decrease in this density with radius. The velocities 
suggest that the C2-emitting gas fills the space between 
the two walls, and the narrower width of C2 seen in the 
centre of the channel implies that the interaction between 
the hot and cool gas is weaker here. For a temperature of 
T = 10^ K (characteristic of the gas associated with the 
hot wind fluid) the sound speed is ~500 kms^^, meaning 
that the expansion of the channel walls is strongly sub- 
sonic (even after taking into account projection effects). 
This is not unexpected since the surrounding medium 
is known to be at high pressure dW073 . Paper II). The 
identification of this outflow channel therefore provides 
us with a rare opportunity to gain a deeper insight into 
the local physical conditions which lead to superwinds, 
and will be discussed more fully in Paper II. 

5.3. The relationship between the narrow and broad line 
components 

We can now examine in more detail what the rela- 
tionship between the narrow and broad line components 
might be, and why C2 appears to rotate at a slower ve- 
locity than CI. 

If C2 represents turbulent gas stirred up by the inter- 
action between the hot, fast-flowing wind and cooler gas, 
we would expect C2 to become more dominant as the 
gas becomes increasingly perturbed and entrained into 
the outward flow. Thus, if this process is occurring over 
the length of a channel such as we describe above, then 
the emission of CI will be weighted towards the inner- 
most region (the base of the outflow channel), whereas 
the C2 emission will be weighted towards regions fur- 
ther out along the flow (i.e. extra-planar regions). The 
velocity gradient difference can therefore be explained, 
since when observing a highly inclined galaxy, gas lo- 
cated further out in the disk would exhibit a shallower 
major axis line-of-sight velocity gradient because here the 
radial (line-of-sight) component of the velocity vector is 
smaller. The precise situation, however, is complicated 
by the presence of the bar in M82, since a bar can give rise 
to a variety of major axis rotation curves (including solid- 
body-like, or steeper) depending o n the bar orientatio n 
with respect to the major axis PA (jD icaire 
The outflow channel model also helps explain why the 
CI radial velocity maps of the inner wind regions show 
a large number of localised features that appear to be 
mirrored by C2 (e.g. see pseudo-slit 'b' in Fig. [T^ since 
the two components are naturally dynamically linked. 
As the tube expands with radius and continues to ro- 
tate with the galaxy disk, the outflow speed of the gas 
will decrease du e to con servation of angula r momentum 
(jSeaauist & Cla rk 200 li I Walter "erni[200l IC^evc 2004). 
If the outflow speed is slow enough then the channel will 
become distorted by these effects; if it is fast enough 
then the effects will not be significant and the channel 



will remain straight. Since the filament in our observa- 
tions appears coherent and linear, we can infer that the 
outflow channels must be reasonably short lived and con- 
tain gas with flow velocities that are large compared to 
the local rotation speed. 

5.4. Implications for the superwind system 

Our CMOS observations sample the gas conditions at 
the base of the superwind in unprecedented detail. Here 
the wind is clearly dominated by processes occurring 
on small spatial scales; these structures are associated 
with the compressed, fragmentary ISM and the SSCs. 
Within this inner region, the flow is not coUimated as is 
seen in the outer wind, but fairly dynamically chaotic. 
Here, fragments of the ISM represent the remnants of 
interstellar clouds as well as shells or bubbles, shredded 
by multiple interactions and shocks. Clearly, however, 
some of this gas is ordered into coherent outflow chan- 
nels, presumably formed throu gh cluster wind interac- 
tions (jTenorio-Tagle et"aI]|2003D as described above. 

Finding C2 emission associated with turbulent clump 
surface layers and outflow channels throughout the disk 
of M82 provides powerful direct evidence for the exis- 
tence of mass- loading, not just at the wind base, but over 
a large, spatially extended area extending far into the in- 
ner wind. The need for large-scale mass-loading in the 
wind of M82 has been inferred ever sin ce the first analytic 
wind models (jChevalier fc Cleg^llQSSl) were compa red to 
X-ray measurements (see also lSuchkov et al"]ll994[ ). 

The near coexistence of multiple phases of the ISM 
should lead to complex cooling processes. For example, 
the measured expansion velocities of the ionized gas fea- 
tures are much less than the the sound speed in the hot 
gas, implying that the mixed material at the base of the 
wind is strongly subsonic. These will range from mod- 
erate cooling as indicated to fit the b ulk X-ray intensity 
of the wind (e.g. ISuchkov et al.|[T994f ) to cases where the 
dense ISM might locally completely dominate and cool 
the hot phase. We also suspect that some fraction of 
the hot ISM could escape unscathed, and processes as- 
sociated with this hot, low density component will be 
difficult to detect in emission. 

It has been known for many years that as the ISM gas, 
stripped and entrained into the wind outflow, reaches 
greater and greater heights, it retains the signature of 
disk rotation and exhibits a concomitant decrease in ro- 
tation speed. This slowing of rotation above the plane 
was first observed in CO and then in the optical, and 
can be (main ly) accounted for by angular momentum 
conservation llSofue et al.l 11991 ISeaquist fc Clarkll2001l : 
IWalter et aLll2002l : lGrevg|2004f ). Our observations con- 
firm this: the imprint of rotation in the velocity patterns 
of CI at all heights is immediately visible in Fig. [13] 
(through the change in sign of velocity relative to sys- 
temic between pseudo-slits 2 and 3 on the east side of 
M82 versus 4 and 5 on the west) and Fig. [181 In agree- 
ment with previous studies, these data also suggest that 
the motions of ISM compressed in expanding bubbles 
along with entrained matter can help to account for the 
larg e mass of extraplanar molecular gas found in M82 
(e.g. iSeaauist fc Claig[200ll : IWalter et al.l[200l . 

6. SUMMARY 
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In this paper we have presented two sets of optical in- 
tegral field spectra of the quintessential starburst galaxy, 
M82. The spatial coverage, depth and resolution of the 
combined datasets have given us an unprecedented op- 
portunity for examining this important galaxy. This pa- 
per is the first of a series examining the optical structure 
of M82's disk and its superwind; here we have focussed 
on the ionized gaseous and stellar dynamics. However, 
the complexity of the dynamics is not easily disentangled 
due to the galaxy's high inclination, strong and variable 
obscuration leading to a degree of scattered emission, the 
bar and starburst-driven outflow. Our main findings and 
conclusions are as follows: 

• The spectral resolution and S/N of our data have al- 
lowed us to examine in detail the shape of the emis- 
sion line profiles observed. In general, they com- 
prise a bright, narrow component (CI; FWHM ^ 30- 
130 kms^^) superimposed on a broad, fainter compo- 
nent (C2; FWHM ~ 150-350 kms'^), with often little 
or no velocity offset between the two. In some regions 
two narrow components are identified where the fainter 
of the two is referred to as C3. Sometimes both nar- 
row components are superimposed on the C2 broad 
emission component. While CI is observed across the 
entirety of our spatial coverage, C2 is identified out to 
a (plane of the sky) radial extent of 1.7 kpc from the 
nucleus, and we do not detect it in most of region B. 

• Within the starburst clumps A and C, both the CI and 
C2 line widths remain very uniform, and we find lo- 
calised ~20-70 kms~^ line splitting in the Ha narrow 
component, which we associate with locally expanding 
shells of compressed, cooler gas that is being pliotoion- 
ized. 

• Narrow CI and C3 component line splitting of 40- 
80 kms~^ continues to larger projected distances. Par- 
ticularly striking is the 50-80 kms""'^ splitting seen 
in region B considering this region is of intermediate 
age and not part of the cur rent starburst (jSmith et al.l 
[2007 : Konstant opoulos et a l. 2008). Possibly the few 
young stars and star clusters in this region suffice to 
produce shell structures, although the possibility of 
other mechanisms should not be excluded. 

• Moving out into the inner wind, we find that the line 
splitting increases to >100 kms^^ and both the line 
splitting and FWHM morphologies begin to follow mi- 
nor axis oriented patterns. This suggests that here 
some of the gas (particularly the C2-emitting mate- 
rial) is associated with linear streaming filaments or 
shocked interfaces. This picture receives further sup- 
port from the discovery of a distinctive outflow channel 
described in Section [5?2l 

• On average, the narrowest C2 lines are found near the 
disk midplane. The widths of the C2 line then increase 
with height above and below the disk to the edges of 
our spatial coverage. This trend is seen in both Ha 
C2 and [S ill] CI, and is an indication that emission in 
the broad component becomes more important with 
increasing height within the energy injection zone. 

• The imprint of disk rotation is seen at all heights mean- 
ing that orbital angular momentum is being carried out 



into the wind (cf. lGrev3l2004f ). The shallower rotation 
amplitude observed for C2, however, is consistent with 
it arising at larger heights where conservation of angu- 
lar momentum has resulted in reduced rotation speeds. 
This pattern has been previously seen in molecular and 
ionized gas emission lines observed in the radio where 
galaxy transparency is not a major issue, and inter- 
preted as axisymmetric drift. Increased importance of 
the broad C2 optical emission component is also ob- 
served moving outward along the outflow channel. The 
general nature of this phenomenon shows that turbu- 
lent gas mixing into outflowing hot gas becomes in- 
creasingly important below the main disk of M82. 

• The turn-over in the rotation curve observed on both 
sides of the galaxy and large C1-C3 line splitting in 
various places implies that, on average, our line-of- 
sight extends to at least half-way through the disk, 
in conflict with the high levels of obscuration usually 
associated with the nuclear regions of M82. 

• The rotation axis of the gas (as traced by Ha CI and 
C2 and [S ill] CI) is offset from the stellar rotation axis 
(as traced by Can and Ha C4 - the absorption com- 
ponent) and the photometric major axis PA by '^12°, 
not only within the nuclear regions but over the whole 
inner 2 kpc of the disk. This provides further evidence 
of its disturbed nature resulting from its interaction 
within the M81 group. 

• Now that we can identify and map each individual 
line component, it is clear that this central region, 
at the base of the large-scale outflow cone, represents 
a chaotic, complex kinematic environment with many 
overlapping expanding structures located at different 
radii, and that the locations at which the narrow com- 
ponent of Ha begins to split coherently are found to 
vary with distance along the major axis. Only at larger 
radii does the bulk-scale outflow kinematics become 
more coherent. 

• Our [S III] observations suggest that we are are probing 
three distinct gas phases: one originating from the rel- 
atively quiescent ionized gas in or near the star clusters 
(Ha CI and C3); one originating from the same place 
but from a slightly hotter gas phase ([Siii] CI); and 
one from the turbulent mixing layers on the surfaces 
of ISM clumps located in moderate-to-low density re- 
gions exterior to the starburst core in the inner wind 
(Ha C2). 

The next paper in the series will focus on the nebular 
characteristics of the ionized gas in the disk and inner 
wind, including a study of the line ratios, excitations 
and gas densities. At larger radii, we know from deep 
Ha imaging and s pectros copy that a structured, cone- lik e 
outflow develops (|McK9 5: SBH98I: lOh^ama et al.l[200l l. 
Developing an overall morphological picture of the M82 
wind, linking the inner and outer wind regions together, 
will be the subject of Paper III. 
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Fig. 1. — HST / ACS F656N image of the central regions of M82, showing the position of our five 7x5 arcsec IFU fields. The brightest 
starburst clumps are labelled and outlined with dashed lines. The 2.2 ^m nucleus IjLester et al.|[T990f ) is marked by a white cross, the major 
axis by a solid line, and the position of the cluster M82-A1 with a white circle. The white lowercase letters indicate the location from which 
the Ha line profiles shown in Fig.|4]were extracted. 
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Fig. 2. — WIYN _R-band image of M82 (inverted greyscale) overlaid with DensePak continuum flux maps (also inverted) fo r all three 
pointings. Some of the bright starburst clumps/clusters are marked, and the position of the 2.2 fim nucleus IILester et al.|[T990l ) is shown 
with a white cross. The boxes show the position of the GMOS IFUs, and the galaxy major axis (PA = 65°) is indicated by a solid line. 
The letters indicate the spaxels from which the Ha+[Nll] lino profiles shown in Fig.[5]wcre extracted (letter 'b' is located at the end of the 
region A arrow, but has been omitted for clarity). The inset shows an i^5T+WIYN composite (Paper III) covering the same area as that 
shown in the main image, and includes narrow-ljand emission from the wind outflow. The outlines of the DensePak fields arc shown for 
comparison. [A colour version of this figure is included in the on-line version.] 
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Fig. 3. — Example DonsePak spectra from fibre 35 in position 1 showing the full wavelength ranges observed using both gratings. 
Emission and absorption lines are labelled; the continuum shape in the red grating arises because the spectra are not flux calibrated. 
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Spectrum No = 333. Spectrum No = 413. Spectrum No - 302. 




Wavelength (.4) l^avelensth (A) Wavelerglh (A) 

Fig. 4. — Example GMOS line profiles and Gaussian fits; (a) typical of the bright region A (note no line splitting in the narrow 
component); (b) typical of position 3 (close line splitting with a clear broad component); (c) from the south-western corner of position 4 
showing clear line splitting (of 125 kms~^) and the presence of a broad underlying component; (d) from the spot in the north of position 
4 exhibiting very broad and blueshifted C2; (e) typical of the brightest region in the centre of region C (illustrating the consistency with 
which we can fit both the Ha and [Nil] profiles); (f) from the south-east of position 5, showing a narrow component with line splitting of 
'^50 kms^^ and a prominent broad underlying component). The units are in an arbitrary but relative intensity scale. The fit residuals are 
plotted under each panel (in units of cr) with dashed guidelines at ±1 (see text for further explanation). The locations from which these 
profiles were extracted are labelled with the corresponding letters on Fig. [T] 
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Spectruin No - 83. Spectrum No = 47. 
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Fig. 5. — Example DcnscPak Ha and [Nll]A6583 line profiles from (a) the inner- wind in the north of position 1, (b) the centre of region 
A, (c) to the east of region A, and (d) the centre of region B (position 3). This sample was chosen to represent a range in line profile 
shapes including broad, underlying emission features, second narrow peaks, and absorption profiles. The observed data are plotted as a 
solid black line, and the summed model fit is plotted as a solid red line. The flat continuum and individual Gaussian profiles arc plotted as 
red dashed lines relative to the y-axis zero level. The units are in an arbitrary but relative intensity scale. The locations from which these 
profiles were extracted are labelled with the corresponding letters on Fig. [2] 
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Fig. 8.— 



Ho FWHM maps for C3. Labels as for Fig. |6] [A colour version of this figure is included in the on-line version.] 
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Fig. 9. — GMOS Ha radial velocity map for CI. A scale bar is given in kms~^ relative to v^ys (200 kms~^). Four major-axis pseudo-slits 
(0'.'4 width) are represented as lettered lines running north-east to south-west: slit a runs along the galaxy major-axis (PA = 65°; highlighted 
in red). Major-axis slits b, c and d arc parallel to slit a but offset by 2, 4 and 6 arcsecs, respectively. Seven minor-axis pseudo-slits are 
also indicated by numbered lines running north-west to south-east: slit 1 runs along the galaxy minor-axis (PA = 160°; also highlighted 
in red); the widths of these minor-axis pseudo-slits are indicated with dashed lines in each case. Minor-axis slits 2, 3, 4a and b, 5 and 6 
are parallel to slit 1 bu t off set by —6, —3.2, +3.5, +6.0, +9.5 and +12 arcsecs, respectively. The resulting position-velocity diagrams are 
plotted in Figs. [T2l and ll3l [A colour version of this figure is included in the on-line version.] 
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Fig. 12. — Ha position-velocity plots for the four majo r-axi s pseud o-slits shown in Figs [9l [TOl and [TTI The gaseous ([Slll] and PIO) and 
stellar (Can) radial velocity measurements of [McKcith et al.l I I1993I) are plotted in the upper-left panel (slit 'a') for comparison. Radial 
velocities are relative to Veys, and the x-axes show the distance from the minor-axis in both arcsecs (bottom) and parsccs (top; assuming 
1" = 17.5 pc). The radial velocity error bars are for the most part smaller than the symbol sizes, except for C3 where the uncertainties lie 
in the range 10-30 kms~^ (see Section |2.3| | but these have been omitted for clarity. 

[A colour version of this figure is included in the on-line version.] 
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Fig. 13. — H« position-velocity plots for the six mino r-axis pseudo-slits shown in Figs l9l [TOl and [TTI The Ha, [Nll]A6583, and [Slll]A9532 
emission line radial velocity measurements of IMcK95l are plotted in the upper-left panel (slit 1) for comparison. Radial velocities are 
relative to Usys, and the y-axes show the distance from the major-axis in both arcsecs (left) and parsecs (right; assuming 1" = 17.5 pc). 
The vertical dashed lines represent the approximate radial velocity of the gas at the point at which the pseudo-slit crosses the major axis, 
thus showing the zero point with respect to the galaxy rotation. Again, the radial velocity error bars are for the most part smaller than 
the symbol sizes, except for C3 where the uncertainties lie in the range 10-30 kms~^ (see Section |2.3| I but these have been omitted for 
clarity. [A colour version of this figure is included in the on-line version.] 
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Fig. 14. — Ho FWHM (corrected for instrumental broadening) vs. radial velocity (relative to Vbjb) for each component in each IFU field. 
The plots are displayed in approximately the same order as placed on the sky. The dashed ellipse in the position 4 plot encloses points 
associated with the 'knot' discussed in the text. For an estimate of the associated uncertainties, see Section 12.31 [A colour version of this 
figure is included in the on-line version.] 
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Fig. 15. — FWHM (corrected for instrumental contribution) maps for the Ha CI and C2 emission components, [S ll] A6717,6731 and 
[Slll]A9531. The major axis (PA = 65°) is indicated wi th a solid line on the Ho C2, [Sll] C2 and [Slll] CI plots. The locations of the two 
minor axis position-velocity inflection points at ±20" IjMcKeith et al.lll995l ) and discussed in the text arc marked with short thick lines. 
[A colour version of this figure is included in the on-line version.] 
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Fig. 16. — Graph showing the FWHM measurements for positions 1 and 2 (plotted against spaxel number), comparing the widths of all 
the identified components of Hq (black), [Sll] (red), [Sin] (green) and P9 (blue). Although the FWHMs of Hq match the equivalent [Sll] 
widths well (for all three components), the widths of [Slll] CI are, on average, broader. [A colour version of this figure is included in the 
on-line version.] 
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Fig. 18. — DensePak radial velocity maps for all four Ho components (emission: CI, C2, C3, and absorption: C4). The scale is in units 
of heliocentric kms"^ relative to Vgyg. The galaxy major and minor axes are marked with solid lines, and the dashed lines indicate the 
width of the pseudo-slits defined in order to make Figs. I20l and l2ll The dashed line in the top-left panel represents the PA of the gaseous 
rotation axis (determined by eye). [A colour version of this figure is included in the on-line version.] 
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Fig. 19. — DensePak radial velocity maps in [Slll]A9531 (CI and C2; upper two panels), P9 (lower-left panel) and Can (lower-right 
panel). The scale is in units of heliocentric kms~^ relative to Vgyg. The galaxy major and minor axes are marked with solid lines, and the 
dashed lines indicate the width of the pseudo-slits defined in order to make Figs 1201 and 1211 [A colour version of this figure is included in 
the on-line version.] 
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Fig. 20. — Major-axis position-velocity diagrams f or th e n ebu la emission lines (left) and stellar absorption lines (right), extracted from 
the corresponding 7" wide pseudo-slit shown in Figs 1181 and 1191 Velocities are heliocentrically corrected and shown with respect to v^ys . 
On the left panel error bars are not shown, since in most c ases they are ap proximately the same size as the symbols. -|- symbols represent 
the optical [Nn]A6 583 mea surements from [Castles et al.l I I1991I . ICMG91T1 and the near-IR gas emission line and stellar absorption line 
measurements from lMcK93l . as indicated. [A colour version of this figure is included in the on-line version.] 
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Fig. 21. — Minor axis position- v eloc ity diag rams for Ho (left panel) and [S m]A9531 and P9 (right panel), extracted from the corresponding 
7" wide pseudo-slit shown in Figs ll8l and ll9l Velocities are heliocentrically corrected and shown with respect to v^yB- Error bars have been 
omitted for clarity; in most cases they are between 5—20 km s^^ with the larger errors applying to the fainter line components. The optical 
and near-IR gas emission line measurements from lMcK95l arc plotted with + symbols. [A colour version of this figure is included in the 
on-line version,] 



